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We discuss the extension of the matrix element method (MEM) to Next-to-Leading Order 
(NLO) in perturbation theory. In particular we focus on the production of a Standard Model 
Higgs boson which decays into four leptons. 
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1 Introduction 

The matrix element metho d u ses fixed order matrix elements to calculate probabilities for ex- 
clusive experimental events^21. By varying the underlying theoretical parameters of the matrix 
element one can determine the best fit values between theory and data. The ensemble data set 
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can then be used to define a likelihood associating the input parameters with the experimental 



data set. By varying the underlying theoretical parameters used in the matrix elements one 
can obtain multiple likelihoods, the maximum likelihood corresponds to the best fit parameters 
linking the theory model to data. The MEM has been extensively used in experimental analyses, 
see e.g. ref.^ for a review of the MEM's application to the measurement of the top mass. This 
t^j- . talk illustrates how this can be extended to NLO in perturbation theory^. 

m 

CO 



2 The Matrix Element Method at LO and NLO 



The primary difficulty inherent in modeling experimental events with fixed order matrix elements 
occurs when attempting to map an experimentally observed set of objects p to a Born phase 
space point p in which the beams are along the z-axis. We define the sum over the particles 
identified with the Born final state as X, i.e. X = —J2?=iPi- For a generic event X x and 
X y 0, which is incompatible with our assumption that the initial state partons are aligned 
with the beam. In order to overcome this obstacle we perform a Lorentz transformation to a 
frame in which X? = 0. This preserves all of the Lorentz invariant quantities associated with 
the experimental event. We now need to construct the longitudinal components of the initial 
state particles which are fixed through the corresponding components of the final state particles. 
However the Lorentz boost which we performed does not uniquely fix these components. In other 
words, there are multiple frames in which the final state particles are balanced in px connected 
to each other by longitudinal boosts. We refer to this collection of frames as the MEM frame. 
In order to provide an unbiased weight we must integrate over all allowed boosts. We note that 
the matrix element is a Lorentz scalar and as such the only boost dependent term we need to 
consider for the MEM is the integration over parton distribution functions 

£ij (Sabi 3Cl i -Eli) — f dx a ■ ■ . (1) 

Jxi SX a S a l) 

Combining this boost integration with the boost invariant matrix element B 13 allows us to 
construct the probability density function for the MEM accurate to LO, 

V{-x\Cl) = ~To I d y J ^ij( s ab,xi,x u )B t ^(p a ,p b ,y)W(yi,y) . (2) 



Here W(x, y) represents the experimental transfer function which models the detector effects. 
We will assume that W(x, y) = <5(x — y), which is valid for identified muons and electrons. 

In order to extend the MEM formalism to NLO we need to incorporate both virtual and 
real contributions into the weight under the constraint that the weight should be evaluated for a 
fixed experimental input event. We imagine that we have performed the Lorentz boost described 
above such that the experimental event has the kinematics of a Born phase space point x. In 
this setup our NLO calculation should be formulated as follows, 

j n NLO( \ 

\ (J =^(x)+V n (x). (3) 

That is, we define the virtual Vh(x) and real -Rq(x) parts of the calculation separately as a 
function of the Born phase space point x. Summing over the Born phase results in the usual 
NLO cross section. Defining the virtual phase space is straightforward since this piece shares 
the same phase space as the Born contribution, the virtual piece is thus, 

Vn(x) = Cij(s ab ,xi,x u )(B^{p a ,p b ,y:) + V^(p a ,p b ,x] 

+ dz(v m (z,x)®£ m (z,s ab ,xi,x u )] B$(p a ,p b ,x.). (4) 

m=0 J ^ ' ij 

Here the first line represents the contributions from the Born matrix element and the virtual- 
born interference terms which occur at one- loop V$f . These pieces contain divergences which are 
regulated through a subtraction scheme. These subtractions are denoted in the second line and 
factor onto the Born matrix element. We observe that since we are considering electro- weak final 
states the subtractions are for initial state singularities. This results in convolution integrals 
between the dipole parameter z and the boost integration. This is shown schematically by the 
sum over m in the above equation. 

In addition to the virtual contributions we must also define the real corrections associated 
with the radiation of an additional parton. These pieces are more troublesome since they reside 
in a higher dimensional phase space than the Born. In order to maintain the desired mapping 
to the Born phase space point we use a Forward Branching Phase Space generator (FBPS)^, 
this provides the following factorisation. 

d$( Pa + p b -^Q + Pr ) = d<S>(p a +p b -> Q) x d$ FBPS (p 

aiPbiPr) * & veto ' (5) 

Hatted momentum represent an underlying Born topology whilst the un-hatted momenta are 
the real phase space point. We note that the observed particles Q are identical to their Born 
counterparts. In terms of the kinematic invariants the FBPS is given by, 

d$ FBPS (p a ,p b ,p r ) = j^-^ (^f) dt ar dt rbd<f> ■ (6) 

Using the FBPS we can now explicitly define Rq(x) as, 

Rn(x) = J d^ FB p B (pa,p b ,p r )(Cij(Sab,Xi,X u )TZ^(p a ,p b ,X,p r ) 

- ]T djisat,, xf, x™)D m (p a , Pb , Pr )B%(p a ,p b , x)) . (7) 



The first term in the above equation represents the integration over the FBPS of the real matrix 
elements TV 3 . In certain regions of phase space these terms develop singularities which are 
regulated by the subtraction terms defined in the second line of the equation. 
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Figure 1: The log- likelihood difference for background only and signal plus background, for a Higgs boson search 
in the channel, H — > ZZ* — > 4 leptons. Positive values of the difference indicate that the background-only 
hypothesis is more likely than the signal plus background one. The blue and magenta lines represent the 1- and 

2-er limits respectively. 



We are now in a position to define the NLO probability density to be used in the MEM, 

P(x|fi) = — (V n (x) + i^(x)) . (8) 

We have suppressed the dependence on the transfer functions, assuming perfectly resolved par- 
ticles. In the next section we will present an application for which this assumption is reasonable, 
namely the production of a SM Higgs boson and its decay to four charged leptons. The future 
applications of the method for LHC physics are widespread. 



3 The Search for the SM Higgs boson 

The search for the Higgs boson is one of the most pressing in experimental particle physics. 
Current LHC limits indicate that, if it exists, then the SM Higgs has a mass in the range 120- 
125 GevEd. One of the most promising decay modes in which to extract the Higgs properties 
is the decay of the Higgs to ZZ which subsequent decays to charged leptons ZZ — > AO,. In this 
instance the final state is fully r econ structed and contains particles which the general purpose 
detectors can measure accurately^El. In this example we generate samples of unweighted events 
produced from a NLO sample, directly in the MEM frame. The underlying physics is identical 
to that implemented in MCFM^l. We assume that no Higgs boson exists and proceed to set 
limits using the MEM. In Fig.[T]we present a results from a single pseudo-experiment for around 
250 events. In this example we sweep over a range of Higgs masses and set limits. NLO sets a 
limit of 100 < m H < 430 GeV, whilst LO sets a limit of 120 < m H < 380 GeV. In Fig. [2] we 
generate multiple pseudo experiments and test a single hypothesis (mn = 200 GeV). As such 
we are able to discern the differences between LO and NLO in a more systematic nature. We 
observe that in general there are observable differences between LO and NLO. The NLO results 
set better limits, however this is hardly surprising given that the underlying sample is NLO. 



4 Conclusions 

We have illustrated how the matrix element method may be extended to NLO in perturbation 
theory. As an example we have considered the decay of the SM Higgs boson to four charged 
leptons. 




A=Log(L B /L s+B ) 



Figure 2: Pseudo experiments testing the hypothesis that there is a Higgs boson with tuh = 200 GeV. We generate 
pseudo-experiments which consist only of background and no Higgs signal. As such the most common outcome 
is that the signal plus background hypothesis is less likely than the background only. 
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